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The determination of the masses of Supersymmetric particles such as the selectron, for energies 
above threshold using the energy end-points method is subject to signal deconvolution difficulties and 
to Standard Model and Supersymmetry backgrounds. The important features of ^ production 
are used to design an experimentally robust method, with good resolution, both for determining 
the cr^l and the x? masses and for suppressing backgrounds. Additional features, such as the 
determination of the relative leptonic branching ratios of the selectron are present in the method. 
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The determination of Supersymmetric particle masses 
using the energy end-point method |Q] is well known. The 
measurement of selectron masses is subject to two exper¬ 
imental difficulties: on the one hand the energy distribu¬ 
tion of the visible particles in the event is an overlap of 4 
“box like” distributions due to the production channels 

other, the Supersym¬ 
metry (SUSY) signal is masked by very large Standard 
Model (SM) backgrounds, such as W+W“ and 7*7*. The 
consequences are, for the former the difficulty of resolv¬ 
ing overlapping “box”-edges in the lower energy region 
making it hard to determine which edge is due to which 
SUSY particle decay, while for the latter the masking of 
the SUSY-signal by large SM backgrounds. 

Many other studies |^|-^ of the determination of su¬ 
persymmetric masses via the energy spectrum of the ob¬ 
served particles have been carried out showing the useful¬ 
ness of the technique and the levels of accuracy possible. 
The complications in the selectron spectrum are being 
solved here. 

During Snowmass-2001 we realised that the difference 
of the observed positron and electron (from the cl.r de¬ 
cay) energy distributions enhanced by the difference in 
cross sections with incident electron polarization (possi¬ 
ble only with a Linear Collider) can solve these problems 
and provide a method by which we can determine the 
masses unambiguously with good resolution. In addition 
it offers new features, built in redundancy, the deter¬ 
mination of the partial leptonic branching ratios to the 
channels xfve and XiC. Information about the xt and 
masses may, in principle, also be determined. 

The e'^, e~ energy distributions difference is given by: 

A{E) = Lumi ■ (auL - ^lh) ' (1) 

where, for instance, cflr is the production cross section 
for eje)j, while R'{E) and L'{E) are the appropriate 
energy box distributions (for the incident electron right- 
handed and left-handed polarization) each normalised to 


unity. The method presented here removes the contribu¬ 
tion to the energy spectrum from the reactions producing 
and eje)), which solves the mass measurement dif¬ 
ficulty and reduces the mass errors. 

Asymmetric boosts are present when and are 
produced in the same reaction. The values of the boosts 
being in this case: 
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The energies of the electron and positron which are the 
particles visible in the event are related to the decay CM 
energy. Since rrie «C we obtain the lower and higher 

bounds of the electron/positron energy distribution in 
the LAB reference frame: 


U,±(iJJ)=A:±(7+v'7'-l) 


(LO) = E:^ (7 - v7^) = + V7^) (3) 

Solving for E*^- and 7 we obtain: 

El = lE,±{HI)-E,±{LO) 
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where A = 7 4- H — 1, and 7 = A/2 -|- 1/(2A). We can 
determine the relationship between the mass errors and 
the energy-edge errors. To a very good approximation 
this follows from: 

_ heL,R _ ^Il,r ~ 1 /<5Ag^,R\ ..s 

M~eR.R 1,R,R + 1 1 Ag,,„ ) ^ ^ 

This becomes to a good approximation 
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where we assume that 5{Ee)/Ef. is the same at the low 
and high end of the energy spectrum. 

The complete set of equations for determining the 
masses from the energy spectrum end-points is: 

= y/s + - M|^) 

= yiV {Ml - Ml) 


2El^M,^ = Ml - Ml 

2El^M,, = Ml - Ml (7) 

The system, although quadratic, allows one to substitute 
the quadratic terms and reach a linear solution: 
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The fit process that searches for the masses Mgj^, Mg^ 
and M^o is carried out with MINUIT [Q that obtains the 
best solution to the set of equations above. 

The properties of the above defined distribution are 
quite simple and robust. First, there are only 2 (well 
separated) “boxes” as seen in Fig. Q one positive, one 
negative, with an overlapping mid-region. 
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FIG. 1. Energy distributions difference between e”*" and e~ 
for only selectron decays. Using the difference of the differ¬ 
ence, between RightgQ% and Leftso% e~ beam polarisations, 
the differences are enhanced and any detector asymmetries 
between charged tracks are effectively cancelled. 


the edge-to-particle assignments are clear, the positive 
part of the distribution corresponding to e^. Thirdly, 
since all SM backgrounds have the same energy distri¬ 
bution for the visible e'^ and e“, their difference does 
not exhibit any contribution from Standard Model back¬ 
grounds. Even if there is a detector asymmetry between 
positive and negative tracks, this is eliminated through 
the use of the polarised version of the above difference as 
described in eqn. 1 and shown in Fig. 1. 

It should be noted also that even though the back¬ 
grounds are reduced to the level of their statistical fluctu¬ 
ations, they can still be large and affect the signal. Hence 
removal of the background to first order by kinematical 
cuts or strategically located detector vetoes is still useful. 
It has been shown that the 7 * 7 * background can be 
suppressed to almost zero using these techniques. These 
cuts fail however when ^ and are close in mass, 
for this case a Very-Forward Veto device is crucial. The 
method proposed here greatly alleviates the kinematic 
severity of the cuts needed to control this background. 


1000 

800 

600 

400 

200 


"0 20 40 60 80 100 120 140 160 180 200 GeV 

FIT = Pj 

FIG. 2. Parametrisation of the secondary leptonic decays of 
the Xi^ ■ The chosen model was an exponential, which requires 
one shape parameter. 

The procedure is complicated slightly by the fact that 
for a region of the SUSY parameter space has two lep¬ 
tonic decay channels, one involving and one XiC^. 

It is assumed that the secondary hadronic decays can 
be eliminated, however the leptonic decays with x^Ve 
have the same signature, an in the final state as the 
preferred decays and have to be accomodated in 

the fit. The difference distribution becomes in this case: 

A{E)= Lumi ■ {(Jrl - (^lr) ■ (/ + /') ■ 
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This means that the edges can be comfortably resolved, 
as none of them overlap in any energy range. Secondly, 
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where X'{E) is the energy distribution of the visible 
from the secondary yf decay, and / and /' are 
the branching fractions to L'{E) and X'{E) respectively. 
The power of the method in its initial form is that this 
histogram has to be normalised to zero, and hence any 
pedestal can be determined and subtracted. It also has 2 
normalisation conditions related to the R'{E) and L'{E) 
boxes, that connect the edge positions with the box- 
heights. To first order there are only 3 free parameters in 
the fit: and M^o. This provides a very robust 

fit. In the present case, the X'{E) distribution (figure |^) 
has to be parametrised, and this adds an extra parameter 
to the fit. The chosen model was an exponential shape 
for X'{E), hence only one extra shape parameter. The 
general idea of the simple fit remains: a robust fit (figure 
^ with few, tightly bound parameters. The fit yields also 
the //(/ + /') and /'/(/ + /') relative leptonic branch¬ 
ing fractions for the xt^e and y^e. The X'{E) shape 
parameter gives general information about the xf and 
masses. 
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This analysis is carried out for the MSSM parame¬ 
ters mo=100 GeV, mi/2=250 GeV, Ao=0, tan(/3) = 10, 
^^positive. The resultant masses (in GeV) of the super- 
symmetric particles in our study is mg^ = 143 , mg^=202, 
m^o=96, m-±=174, and m~±=135. The errors in the 

Xi 'i 

mass determination versus luminosity are shown in fig¬ 
ure H for the “idealistic” case with no x^t'e channel, as 
well as for the “realistic” case including the xt^e- An in¬ 
dication of saturation of the mass-errors with luminosity 
is shown by the fit. The luminosity at which the error 
reaches 80% of its limiting value gives a range of 500 - 
1000 fb~^ as the saturation limit. It can be seen that the 
largest difference in mass-errors between the “ideal” and 
“real” case is for due to its corresponding energy- 
“box” neighboring the perturbing X'{E) distribution. 
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FIG. 4. Mass determination errors versus luminosity, for 
the prototype distribution that includes only the Xie decay 
and for the realistic one including also the xf t'e decay. An 
indication of saturation with luminosity is seen. 


FIG. 3. Fit to the electron energy distribution includ¬ 
ing the xf decay mode of the El based on an in¬ 
tegrated luminosity of 2000 fb“^ and a collision en¬ 
ergy of 500 GeV. The input masses (in GeV) are 
Msj, = 142.96, = 202.07, M^o = 95.70. The fit gives 

= 143.12±0.23, Mg^ = 202.27±0.11, M^o = 95.72±0.12. 


The method is very robust and stable, its results be¬ 
ing not only the masses involved, but also possibly the 
relative leptonic branching ratios, and some tangential 
information about the xf and rf' masses. This last as¬ 
pect is currently under study. The resolution obtained is 
comparable to that obtained with production threshold 
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studies A study of the mass determination errors 

versus luminosity indicates a possible saturation of the 
errors from the I/VN law, occuring when the luminosity 
is of the order of 500 - 1000 fb~^. 
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